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INTRODUCTION 


This  paper  discusses  the  Phillips  Laboratory's  Electric 
Propulsion  Laboratory  (PLEPL)  pulsed  magnetoplasmadynamic  (MPD) 
thruster  performance  experiment.  Goals  in  measurement,  model 
selection,  calibration  procedure  and  sources  of  error  are  de¬ 
scribed.  Also  included  is  a  short  discussion  of  statistical 
error  analysis  of  the  data  collected. 

Background 

Certain  performance  characteristics  are  desired  for  compari¬ 
son  and  evaluation  when  demonstrating  new  technology  rockets. 
Thrust  and  specific  impulse  are  among  these  characteristics. 
Similarly,  electric  power  use  efficiency  is  used  to  evaluate 
electric  thrusters. 

An  MPD  thruster  uses  electromagnetic  fields  to  accelerate  an 
ionized  gas  (plasma) .  Present  experimental  MPDs  have  specific 
impulses  of  2000  -  4000  seconds  and  efficiencies  in  the  25  -  30 
percent  range  (Refs.  1,2).  High  thrust  density  and  structural 
simplicity  make  MPD  thrusters  attractive  for  a  variety  of  future 
space  missions,  such  as  spacecraft  maneuvering  and  orbit  transfer 
(Ref.  3). 

There  currently  exists  little  performance  baseline  data  on 
long  electrode  length,  multi-megawatt  pulsed  MPD  thrusters  (Refs. 
4,5).  The  PLEPL  is  presently  testing  a  5  cm  diameter  anode  (half 
scale)  variable  geometry  MPD  thruster.  PLEPL 's  goal  is  to  map 
pulsed  MPD  thruster  performance  as  a  function  of  electrode  geome¬ 
try.  The  end  goal  is  to  provide  data,  curves,  and  trends  for 
different  anode  and  cathode  length  configurations.  Semi-empirical 
theoretical  models  may  then  be  derived  from  the  data. 

For  pulsed  MPD  thrusters,  a  high  level  of  electrical  noise 
limits  the  use  of  traditional  thrust  measurement  techniques  (Ref. 
6) .  Thrust  from  pulsed  MPD  effect  lasts  a  few  milliseconds  and 
is  on  the  order  of  tens  of  pounds.  Direct  thrust  measurement  is 
difficult  because  of  the  high  currents  and  voltages  associated 
with  electric  propulsion. 

The  PLEPL  built  a  total  impulse  stand  to  measure  impulse 
(rather  than  thrust) .  Pulsed  MPD  thrust  is  found  by  dividing  the 
impulse  by  the  firing  time.  The  stand  consists  of  a  dial  indica¬ 
tor  attached  to  the  thruster  with  the  measurement  end  placed 
against  a  rigid  support. 


Project  Goals 

The  goal  of  this  project  is  to  map  MPD  thruster  performance 
as  a  function  of  electrode  geometry.  Theoretical  models  and 
experiments  have  revealed  a  dependence  of  thrust  on  the  ratio  of 
anode  diameter  to  cathode  diameter  (Ref.  7) .  There  has  also  been 
work  on  scaling  of  thruster  performance  (Ref.  8) .  Finally,  some 
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geometry  performance  mapping  has  been  conducted  by  Ref.  9.  This 
project  intends  to  extend  the  mappings  to  variation  in  cathode 
and  anode  length.  Thus,  the  cathode  and  anode  length  influence  on 
MPD  performance  may  be  determined,  resulting  in  a  better  under¬ 
standing  of  MPD  physics  and  a  knowledge  of  the  best  performing 
geometry . 


FACILITY  COMPOMEMTS 

The  pulsed  multi-megawatt  MPD  facility,  also  known  as  Cham¬ 
ber  #2,  consists  of  an  8  foot  diameter  by  12  foot  long  stainless 
steel  cylinder.  The  chamber  has  one  door  and  six  view  ports  along 
the  sides  (Figs.  1,2} . 


Figure  1 
Facility  Layout 


Fgure  2 

Pulsed  MPD  Thruster  Test  Configuration 
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A  Stokes  mechanical  pump  (Model  412H-10) ,  a  Roots  blower 
(Model  615  RGS) ,  and  two  Varian  12  inch  diffusion  pumps  (Model 
0185)  bring  chamber  pressure  down  to  5E-5  Torr  with  no  propellant 
flow  (Fig.l).  The  pumps  also  allow  a  return  to  this  pressure 
within  5  minutes  after  a  0.1  second  gas  pulse  and  a  6  grams  per 
second  mass  flow.  Heavier  gases  such  as  Xenon  may  take  as  much  as 
10  minutes  to  recover  to  the  minimum  pressure. 

The  propellant  system  consists  of  a  326  cubic  foot  argon  T- 
bottle  and  regulator,  with  28  cubic  foot  bottles  of  Hydrogen, 
Helium,  Neon,  Xenon,  and  Krypton  also  plumbed  into  the  system.  A 
thermocouple  attached  to  the  propellant  line  outside  of  the 
chamber  gives  a  measure  of  the  propellant  temperature.  The  pro¬ 
pellant  lines  feed  into  a  231  cubic  inch  plenum  located  inside 
the  chamber.  At  the  plenum,  the  mass  flow  is  calculated  by  meas¬ 
uring  the  pressure  drop  per  gas  pulse. 

A  Del  Electronics  Corp  power  supply  (Model  HPS-1-8000-3 ) 
supplies  up  to  8  kWe  of  power  to  a  pulse  forming  network  (PEN) . 
The  PFN  is  a  ten  section  LC  network  with  a  nominal  O.Ol  ohm 
output  impedance.  Each  section  consists  of  three  2000  microfarad 
Maxwell  capacitors  (Model  33800)  connected  in  parallel  and  a  5 
turn,  0.53  microhenry  inductor.  Together  they  release  a  one 
millisecond  current  pulse  at  up  to  40  kAmps  and  400  volts  to  the 
thruster,  assuming  a  perfectly  matched  load  (Figs.  3,4). 
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Rgure  3 

Schematic-Pulse  Forming  Network 
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Rgure  4 

Pulse  Forming  Network 
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Facility  data  is  collected  with  a  variety  of  devices.  Cham¬ 
ber  pressure  is  measured  with  a  Varian  843  vacuum  ionization 
gauge.  Pearson  Electronics  current  pulse  transformers  (Model 
301X)  measure  the  current  flowing  through  the  thruster.  Voltage 
across  the  thruster  is  measured  by  probes  attached  to  the  PFN. 
Both  voltage  and  current  data  are  collected  on  a  Tektronix  DAS 
601  digital  oscilloscope. 


MSASUREMEMT  MODEL 

The  thruster's  motion  is  modeled  for  experimental  simplicity 
and  repeatability.  The  thruster  is  a  copper  cylinder  (about  116 
pounds)  suspended  on  two  nylon  ropes  (Fig.  5) .  Operating  a 
steady-state  MPD  thruster  and  knowing  the  exact  thruster  weight 
and  displacement  can  yield  a  very  accurate  thrust  measurement. 
However,  the  PLEPL  thruster  operates  in  a  pulsed  mode  with  a 
total  firing  time  of  1.5-2. 5  milliseconds.  This  negates  any 
attempt  to  use  a  static  force  vector  model.  The  alternative  is 
to  model  the  thruster  as  a  pendulum.  An  instantaneous  momentum 
change  is  a  good  approximation  (Eq.  1) ,  since  the  firing  pulse 
time  is  much  less  than  the  period  of  the  thruster  swing. 


Impulse  »  mass  *  change  in  velocity  Equation  1 


Figure  5 
MPD  Thruster 


The  calibration  model  selected  is  an  inelastic  collision  of 
a  pendulum  colliding  with  another  pendulum.  (Fig.  6) (Ref.  10) 

Let  m  =  mass  ot  ball 

M  =  moss  of  thruster 

V  =  velocity  of  boll 
b 

=  velocity  of  thruster 

i<  =  coefficient  of  restitution 
(value  between  0  and  l) 

subscript  1  indicates  values 
before  collision  and  subscript  2 
values  after  collision. 
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L  =  pendulum  lenmji 
s  =:  honzorrtal  ball  displacement 
h  =  vertical  ball  displacement 
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Rgure  6a 

Ball  Impacting  Stand 


support  ropes 
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Rgure  6b 
MPD  Pendulum 


From  conservation  of  momentum  (inelastic  collision),  the 
equation  is: 


k(mV  +  MV  )  =  rnV  +  MV 
^  bl  tr  c2  t2 


Since  the  velocity  of  the  thruster  is  zero  before  collision. 


/ 

ti 
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Equation  2 


5 


This  reduces  the  conservation  of  momentum  equation  to: 


=  nnV^2+ 


Equation  3 


Solving  Equation  3  for  velocity  of  ball  after  collision  gives: 


kmVbi  -  MVj,2 


b2 


Equation  4 


m 


From  conservation  of  energy,  we  have  the  equation: 

2  .  2 
4-  ..  . 

2’  t2 


Substituting  in  =  0  from  Equation  2,  this  reduces  to: 


kmV  =  mV  +  MV 
b1  b2  t2 


Equation  5 


Substituting  Equation  4  into  Equation  5  and  solving  for  the 
velocity  of  the  thruster  after  collision  yields: 


Equation  6 


lc=*l  (elastic  collision),  Equation  6  reduces  to: 
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The  velocity  of  the  ball  before  impact  can  be  determined  by 
equating  its  potential  energy  before  release  to  its  kinetic 
energy  at  the  bottom  of  its  swing. 

1  2 

mqh  =  TmV 
‘  bi 

Thus  the  total  impulse  imparted  to  the  thruster  is: 
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where  g  is  acceleration  due  to  gravity,  and  h  is  the  vertical 
distance  the  ball  traveled  from  release  to  the  bottom  of  its 
swing.  Using  trigonometric  ratios  (Fig.  6a),  the  height,  h,  can 
be  equated  to  the  length  of  the  pendulum,  L,  and  the  horizontal 
displacement ,  s . 


h 


L  - 


Equation  10 
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Both  L  and  s  are  much  larger  and  much  easier  to  measure  to  a 
greater  degree  of  accuracy  than  h.  Since  g  is  known  and  m  and  L 
are  measured  for  a  given  pendulum,  it  is  only  necessary  to  meas¬ 
ure  s  to  determine  the  impulse  exchanged  with  the  thruster. 


PHYSICAL  COMPONSHTS 


Total  Impulse  stand 

The  primary  impulse  stand  components  are  a  machinist's  dial 
indicator  attached  to  the  thruster,  and  a  fixed  metal  support  for 
the  dial  indicator  (Fig.  7) .  The  PLEPL  used  two  dial  indicators. 
The  first  one  is  a  Chicago  Dial  Indicator  Co.  indicator  measuring 
in  thousandths  of  an  inch  (part  number  2-ClOO-lOOO) .  It  per¬ 
formed  well,  although  there  are  uncertainties  in  how  vacuum 
conditions  affect  the  indicator  spring  lubricant  and  measurement 
rod.  The  second  dial  indicator  is  a  Checker  brand  indicator  and 
has  a  much  smoother  operation  in  air  and  in  vacuum. 


Figure  7 
Impluse  Stand 
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The  fixed  metal  support  is  a  plate  attached  to  the  vacuum 
chamber  floor  with  a  1.0  inch  Plexiglas  insulation.  The  Plexiglas 
allows  for  smooth  motion  along  the  surface  face.  Repeated  meas¬ 
urements  with  the  Plexiglas  covered  support  have  proven  consist¬ 
ent. 

Video  Recorder  System 

A  video  camera  is  located  outside  the  vacuum  chamber.  The 
video  camera  is  aimed  through  a  view  port  window  at  approximately 
a  45  degree  angle  relative  to  the  dial  face.  During  testing,  the 
dial  deflections  are  viewed  on  a  video  monitor  in  the  PLEPL  control 
center  and  recorded  on  a  video  tape.  The  camera  records  at  a 
speed  of  1000  frames  per  second,  which  is  sufficient  to  find  the 
maximum  deflection. 

Calibration  Component s 

The  calibration  components  consist  of  a  frame  and  attached 
solenoid,  a  rod  for  hanging  calibration  masses,  and  the  masses 
(balls)  themselves.  The  frame  is  an  aluminum  rectangular  frame 
mounted  on  a  flat  base  (Fig.  8) .  A  solenoid  device  is  attached 
to  one  leg  of  the  frame  to  hold  and  release  the  balls.  Also 
attached  to  the  frame  leg,  behind  the  solenoid,  is  a  measuring 
rule  plate  covered  with  graph  paper.  The  rule  plate  serves  as  a 
positioning  index  to  determine  the  initial  horizontal  displace¬ 
ment,  s,  of  the  ball  from  the  thruster  face. 


Figure  8 

Calibration  Components 
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The  balls  hang  from  an  aluminum  pipe  rod  (approximately  0.25 
inch  diameter) .  The  rod  is  attached  to  the  chamber  ceiling  and 
adjusted  such  that  the  balls  barely  touch  the  edge  of  the 
thruster  face. 

The  final  components  in  the  calibration  set  are  the  pendulum 
ball  masses.  Three  balls  are  used,  with  masses  selected  so  that 
the  impulse  imparted  to  the  thruster  is  on  the  same  order  as  the 
impulse  from  a  typical  MPD  test  firing.  Eyehole  screws  are  set 
into  the  balls  and  a  piano  wire  is  attached.  This  allows  a  range 
below  and  above  the  expected  MPD  impulses  for  calibration  pur¬ 
poses.  The  following  table  gives  the  values  of  the  ball  masses 
and  wire  lengths  used. 


Table  1 .  Ball  Parameters 


Bail 

Mass  (grams) 

Length  (meters) 

#1 

46.5  +/-  0.5 

1.1748  +/-0.01 

#2 

169.0+/- 1.0 

1.1684  +/-0.01 

#3 

1  225.0  +/-  1 .0 

1.1862  +/-0.01 

Although  the  masses  remain  constant,  there  are  some  variations  in 
wire  length.  This  happens  if  wires  break  and  have  to  be  replaced 
or  repaired,  therfore  the  wire  length  must  be  remeasured.  The 
length  is  measured  from  the  top  of  the  wire  loop  to  the  center  of 
mass  of  the  ball. 


CALIBRATION  PROCEDURE 

Calibration  is  performed  in  an  open  (vented)  vacuum  chamber 
before  and  after  a  set  of  firings  for  a  given  thruster  geometry. 
This  permits  checking  variations  in  impulse  measurements  and 
gives  an  indication  of  impulse  data  uncertainty  due  to  minor 
changes  in  experimental  configuration. 

Impulse  Stand  Calibration 

The  impulse  stand  calibration  is  a  two  person  operation, 
with  one  person  setting  up  in  the  chamber  and  another  viewing  the 
monitor,  playing  back  the  videotape,  and  recording  the  data. 


Setup  Procedures.  Three  different  horizontal  displacements, 
ranging  from  10  centimeters  to  50  centimeters,  are  used  for  each 
ball.  The  PLEPL  has  used  15  cm,  25  cm  and  35  cm,  roughly,  in 
calibrations.  For  each  ball,  five  identical  drops  are  done  at 
each  displacement  to  ensure  repeatability  and  statistical  analy¬ 
sis.  Forty-five  data  points  are  obtained  in  a  typical  calibra- 
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tion. 


With  the  thruster  in  place,  the  calibration  frame  is  placed 
in  front  of  the  thruster.  The  rod  is  adjusted  so  that  the  ball 
edge  just  touches  the  thruster  face  when  the  ball  hangs  motion¬ 
less.  The  ball  is  also  placed  so  that  it  strikes  the  flat  part 
of  the  face.  The  measuring  rule  plate  is  adjusted  to  roughly  the 
horizontal  distance  desired  for  testing.  The  exact  distance  does 
not  matter,  since  the  ball  distance  from  the  thruster  is  deter¬ 
mined  with  the  measuring  rule.  The  triangular  rule  edge  is 
brought  within  a  centimeter  from  the  thruster  face.  The  distance 
is  measured  and  added  to  the  rule  readings.  When  adjusting  the 
measuring  rule,  care  is  taken  to  make  sure  it  is  parallel  to  the 
chamber  floor  and  to  the  centerline  of  the  thruster. 

Finally,  the  solenoid  is  positioned  to  hold  the  ball,  so 
that  the  ball  swings  freely  and  strikes  the  face  of  the  thruster. 


Calibration  Procedure.  Calibration  begins  with  all  elements  in 
place.  Four  values  are  written  down  at  each  data  point  during 
stand  calibration:  the  initial  dial  reading;  the  maximum  dis¬ 
placement  dial  reading;  the  ball's  initial  horizontal  displace¬ 
ment;  and  the  ball  being  used.  Horizontal  displacements  are  read 
by  looking  at  the  edge  of  the  ball  against  the  measuring  rule 
plate. 

Once  initial  dial  readings  and  horizontal  displacements  are 
recorded,  the  solenoid  is  activated  and  the  ball  released.  Wh.^n 
balls  strike  the  face,  they  return  to  almost  the  same  point  from 
release  (within  2  to  5  centimeters  out  of  15-35  centimeters  or 
usually  about  13  per  cent) .  This  gives  a  measure  of  k,  the 
coefficient  of  restitution.  After  recording  a  portion  of  the 
collision  (15  seconds  after  collision),  the  tape  is  stopped.  The 
tape  is  advanced  frame  by  frame  to  locate  the  maximum  dial  dis¬ 
placement.  The  best  value  of  the  dial  displacement  is  ascertained 
together  with  an  estimate  of  the  uncertainty  in  the  reading. 

The  calibration  frame  and  the  last  ball  are  removed  at 
calibration  completion.  The  chamber  door  is  shut  and  the  vacuum 
chamber  is  ready  for  pump  down.  Once  MPD  thruster  performance 
testing  is  over  and  the  chamber  is  vented  to  atmosphere,  the 
frame  is  brought  back  and  the  same  calibration  procedure  is 
repeated . 


Calibration  Data  Roduction  and  Analysis.  The  calibration  data  is 
reduced  as  follows:  dial  indicator  readings  are  subtracted  f  r 
absolute  dial  displacement;  the  corresponding  total  impulse 
imparted  to  the  thruster  for  those  displacements  is  calculated; 
the  impulse  due  to  the  collision  is  found  by  using  Equations  9 
and  10;  and  a  curve  is  found  correlating  impulse  to  dial  dis¬ 
placements  (with  uncertainty  included) .  Note  that  the  particular 
correlation  is  valid  for  one  set  of  thruster  tests,  i.e.,  one 
geometry.  It  is  expected  that  impulse  as  a  function  of  horizontal 
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displacement  is  a  square  root  function,  as  shown  in  Figure  9. 
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Impulse  Calibration  Curve 
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PULSED  MPO  TESTING 

This  section  describes  those  key  items  involved  in  data 
collection  pertaining  to  thrust  and  efficiency  determination,  and 
the  procedure  used  in  the  PLEPL  pulsed  MPD  in-house  tests. 

Thruster  Operations 

The  half  scale  thruster  is  a  copper  cylinder  (the  anode) 
with  a  2%  thoriated  tungsten  cathode  rod.  A  propellant  tube  and 
power  cables  attach  at  one  end.  A  tungsten  spark  trigger  is 
embedded  in  a  boron  nitride  insulating  plate  inside  the  thruster. 
A  large  capacitor  bank  forms  the  pulse  forming  network  (PFN) .  The 
PFN  creates  a  square  current  pulse  that  is  sent  through  the 
thruster.  When  the  control  center  sends  a  fire  signal,  a  valve 
opens  releasing  a  gas  propellant  pulse.  At  the  same  time,  the 
spark  trigger  releases  a  high  voltage  spark,  which  ionizes  the 
gas  in  the  thruster.  Simultaneously,  the  PFN  fires  a  high  current 
electric  pulse.  This  pulse  creates  the  electric  and  magnetic 
fields  in  the  thruster  that  interact  to  propel  the  ionized  gas 
out  of  the  thruster.  The  resulting  thrust  causes  a  dial  indicator 
deflection  which  is  recorded  on  video  tape.  At  the  same  time,  the 
thruster  voltage  and  current  are  recorded  on  a  digital  oscillo¬ 
scope  (Fig.  10) . 
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Figure  10 

Thruster  Firing  Sequence 


Cold  Gas  Tests 

Before  and  after  MPD  thruster  test  firings,  a  series  of  cold 
gas  tests  are  run  to  determine  the  cold  gas  impulse.  In  theory, 
the  mass  flow  released  through  the  gas  valve  should  be  constant. 
However,  if  the  propellant  line  is  twisted  or  the  propellant  feed 
holes  into  the  thruster  are  stopped  up,  then  there  may  be  a 
change  in  mass  flow.  The  cold  gas  tests  indicate  the  variations 
In  mass  flow,  which  varies  with  different  plenum  pressures.  As 
the  plenvun  is  filled  with  more  gas,  or  as  the  temperature  in  the 
plenvuD  rises,  the  pressure  increases.  This  changes  the  mass  flow 
rate.  Finally,  the  PLEPL  is  interested  in  the  MPD  effects  with 
respect  to  the  impulse.  Since  the  gas  flow  by  itself  causes 
impulse,  the  impulse  due  to  MPD  effects  is  found  by  subtracting 
the  cold  gas  impulse  from  the  measured  impulse  of  an  MPD  firing. 

The  propellant  system  consists  of  a  bottle  of  argon  with 
tubing  leading  to  a  plenum  located  inside  the  chamber.  Transduc¬ 
ers  measure  plenum  pressure  and  temperature.  A  valve  driver 
controls  the  release  of  the  gas  pulse.  Typical  plenum  pres¬ 
sures  are  20  to  30  psia  .  The  temperature  is  read  from  a  thermo¬ 
couple  gauge  connected  to  the  propellant  line  at  the  argon  bot¬ 
tle. 


A  Kistler  force  transducer  (model  202A5/666M1)  measures  how 
long  the  plenum  valve  is  open  when  gas  is  released.  This  data  is 
collected  on  a  Tektronix  digital  oscilloscope  (model  11403) .  The 
curve  generated  is  very  nearly  a  square  wave,  with  a  tail  drop¬ 
ping  off  exponentially.  A  method  to  approximate  the  pulse  time  is 
to  integrate  graphically  and  divide  by  the  height  of  the  nearly 
flat  part  of  the  curve,  thus  getting  the  pulse  width  of  a  perfect 
square  wave  with  an  equivalent  area  under  the  curve.  Although  a 
numerical  quadrature  of  the  curve  data  is  possible,  it  is  unnec¬ 
essary  (Fig.  11) .  A  very  simple  graphical  integration  is  to  take 
the  point  of  the  drop  (visually)  where  the  area  under  the  curve 
past  the  point  is  about  equal  to  the  area  above  the  curve  behind 
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the  point.  This  area  behind  the  point,  if  filled  in,  would  com¬ 
plete  the  square  wave.  Typically,  this  is  about  the  midpoint  of 
the  tail.  Once  this  point  is  found,  the  pulse  time  can  be  read 
directly  from  the  abscissa.  The  PLEPL  has  found  that  the  gas 
pulse  time  is  consistent,  usually  about  106  milliseconds. 


Figure  1 1 

Graphical  Integraton 


The  mass  flow  is  calculated  using  the  ideal  gas  law,  PV=nRT, 
which  is  accurate  for  the  low  pressure,  low  density  argon  used  in 
the  PLEPL  tests.  The  plenum  is  filled  to  the  desired  pressure, 
usually  about  30  psia  (sufficient  time  is  allowed  for  the  gauge 
readings  to  settle  to  a  constant  value) .  With  the  known  volume, 
the  initial  plenum  pressure  and  temperature  are  read  from  gauges. 
The  cold  gas  pulse  is  fired  and  the  pressure  reading  drops.  Due 
to  the  mechanical  inertia  of  the  system,  about  ten  seconds  are 
allowed  for  the  reading  to  settle  to  a  constant  value.  The  gauge 
is  not  allowed  to  settle  too  long,  however,  since  the  pressure 
goes  up  with  temperature,  which  tends  to  rise  over  the  term  of  an 
experiment  (typically  1  degree  Fahrenheit  per  half  hour,  during  a 
four  to  five  hour  test  period) .  From  the  final  pressure  record¬ 
ed,  the  change  in  pressure  is  found  and  change  in  mass  can  be 
correlated.  Dividing  by  the  gas  pulse  time  gives  a  quasi-steady 
state  mass  flow.  The  constants  in  the  PLEPL  mass  flow  calibra¬ 
tions  are  listed  in  Table  2. 

Table  2.  Mass  Flow  Cal  Constants  and  Conversion  Factors 


Experimental  Constants 

Plenum  Volume 

235.3  +/-  0.1  in^ 

Gas  Constant  R 

8.3144  J/mole*K 

Mole  Wtof  Ar 

39.94  g/mole 

Conversion  Factors  i 

Inch^ 

1.6387064E-5  M  1 

Psia 

6.894757E3  N/m 
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The  mass  flow  correlates  linearly  to  initial  plenum  pressure 
(20  psi  is  about  1  g/sec  and  30  psi  is  about  2  g/sec) .  The  MPD 
tests  are  typically  run  at  about  1.5  g/sec. 

Concurrent  with  mass  flow  calibrations  are  cold  gas  impulse 
measurements,  taken  by  recording  dial  indicator  readings  for  each 
cold  gas  pulse.  Typically,  the  displacements  have  been  about  6 
to  8  divisions  of  the  dial  indicator,  although  this  is  dependent 
on  the  alignment  of  the  thruster  and  dial  indicator.  More  impor¬ 
tantly,  the  cold  gas  deflections  have  been  approximately  25%  of 
the  total  deflections  of  MPD  thruster  firings  at  high  (600  V) 
voltages. 

MPD  Thrust  Dstsmination 

From  the  impulse  calibrations,  the  impulse  from  a  firing  can 
be  correlated  to  total  dial  displacement.  This  impulse  repre¬ 
sents  total  impulse  due  to  gas  thrust  and  MPD  thrust. 


Thrust  Calculations.  Since  the  main  interest  is  in  thrust  due  to 
MPD  effects,  the  impulse  due  to  the  gas  thrust  must  be  subtract¬ 
ed. 


=  1  -t-  I  Equation  1 1 

tot  gas  MPD 


The  thrust  due  to  MPD  effects  is  the  MPD  impulse  divided  by 
the  current  pulse  time  taken  from  the  oscilloscope.  This  is 
because  the  electromagnetic  fields  last  only  as  long  as  the 
current  is  running. 


tot 


gas 
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Where  T  is  the  thrust  due  to  MPD  effects  (in  Newtons)  and  t 
is  the  current  pulse  time  (in  seconds).  Specific  impulse  and 
electric  power  efficiency  can  be  calculated  using  the  following 
formulas  (Ref.  11) . 
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Here  V  is  the  voltage  and  J  the  current  flowing  through  the 
thruster.  The  mass  flow  is  m,  and  g  is  the  standard  acceleration 
due  to  gravity  at  sea  level  (9.8066  in/s‘2)  .  The  efficiency  gives 
a  measure  of  how  much  of  the  electric  power  is  transformed  into 
thrust  power. 


Sample  Calculation  of  Thrust.  A  sample  of  the  impulse  data  taken 
so  far  gives  an  indication  of  how  much  a  ••direct”  impulse  meas¬ 
urement  varies  from  the  inelastic  collision  model  impulse.  The 
true  impulse  imparted  to  the  thruster  is  difficult  to  measure, 
but  as  an  approximation,  the  dial  indicator  displacement  can  be 
used  as  an  order  of  magnitude  measure  of  the  thruster  displace¬ 
ment.  This,  along  with  the  calculated  value  of  the  thruster  mass 
(from  volume  of  the  thruster  and  copper  density)  allow  for  an 
approximate  measure  of  the  thruster  impulse.  Table  3  shows  the 
seunple  values  used  and  the  resulting  impulses.  As  can  be  expect¬ 
ed,  the  "direct"  measurement  is  lower  than  the  actual  impulse 
imparted,  which  is  due  to  friction  and  interference  losses  and  is 
discussed  in  the  next  section. 


Table  3.  Example  Thruster  Data 


Bail  Data 

Thruster  Data 

Mass  0.169  kg 

Mass  58.127  kg 

Wire  Length  1 1 9.4  cm 

Rope  Length  1 20  cm 

s  (cm)  H  (cm) 

Dial  DIsp  (cm)  H  (cm) 

14.3  0.859 

28.5  3.45 

38.0  6.21 

0.136  7.7E-5 

1  (N-s) 

1  (N-s) 

0.115 

0.071 

0.219 

0.167 

0.280 

0.226 
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SOURCES  OP  ERRORS 


There  are  a  number  of  errors  and  uncertainties  associated 
with  the  pendulum  impulse  stand,  mass  flow  measurements  and 
electrical  measurements. 

Thruster/  Ball  Support  Errors 

Losses  in  energy  occur  due  to  friction  in  the  calibration 
mass  wires  and  the  thruster  ropes.  A  pendulum  model  assumes  that 
the  ropes  are  identical  and  aligned  in  parallel.  Minor  off-axis 
alignments  cause  the  horizontal  displacement  to  be  less  than  the 
actual  displacement.  The  off-axis  alignment  can  cause  extra, 
complex  modes  of  thruster  oscillations  that  subtract  energy  from 
translational  motion.  Loss  in  energy  due  to  friction  in  calibra¬ 
tion  wires  was  found  to  be  less  than  0.5  percent  in  check-out 
tests . 

The  errors  due  to  changes  in  system  alignment  are  difficult 
to  determine.  However,  the  pre  and  post  impulse  stand  calibra¬ 
tion  give  some  measure  of  the  variation  in  the  calibration  sys¬ 
tem.  As  an  example,  during  calibration  system  check-out,  a  rough¬ 
ly  linear  calibration  curve  was  generated.  The  slope  of  the  line 
was  0.00792  N-s/mil.  After  the  stand  was  removed  and  placed  back 
(more  severe  distortion  than  that  caused  by  vacuum  conditions  or 
firings) ,  calibration  was  repeated.  The  new  slope  measured  was 
0.00688  N-s/mil  (13  percent  change). 

Dial  Zadicatioa  Errors 

Misalignment  of  the  dial  indicator  can  cause  a  lower  dial 
displacement  reading  than  would  normally  be  read.  Ideally,  the 
centerline  of  the  thruster  is  normal  to  the  support  surface.  To 
measure  displacement  along  the  centerline  of  the  thruster,  the 
rod  should  also  be  parallel  to  the  thruster  and  normal  to  the 
support  surface.  Since  both  are  not  always  possible,  the  main 
effort  is  in  placing  the  rod  normal  to  the  support.  This  is 
justifiable  for  correlating  impulse  to  indicated  dial  displace¬ 
ment  since  the  actual  distance  traveled  by  the  thruster  is  not 
easily  measured. 

Video  Tape  Errors 

The  largest  errors  affecting  the  accuracy  of  the  impulse 
measurements  are  those  associated  with  reading  the  dial  indicator 
values  from  the  recorded  video  tape.  The  tape  records  motion  at 
certain  time  intervals.  The  probability  is  small  that  the  true 
maximxim  is  recorded.  The  PLEPL  takes  five  measurements  with  the 
same  ball  in  the  same  position,  so  that  the  errors  in  finding  the 
maximum  are  averaged  out.  Tests  taken  with  larger  numbers  of 
data  indicate  that  five  measurements  usually  result  in  average 
errors  of  less  than  l  division. 

Another  problem  associated  with  the  video  tape  is  that  the 
pictures  vary  in  clarity.  Some  are  quite  clear  and  are  probably 
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indicative  that  a  value  very  close  to  the  true  maximum  has  been 
obtained.  These  can  be  read  to  an  accuracy  of  +/-  0.1  divisions. 
Others  are  blurred  action  shots  that  result  in  uncertainties 
ranging  from  +/-  0.2  divisions  to  +/-  1.5  divisions.  Lighting  and 
shadows  also  make  readings  difficult. 

Other  stand  Errors 

Other  possible  effects  on  the  calibration  are  magnetic  force 
perturbations  and  aerodynamic  drag  on  the  balls.  Although  the 
potential  for  magnetic  interference  is  great  given  the  strong 
magnetic  fields  generated,  tests  indicate  that  they  have  had 
negligible  effect  on  the  impulse  stand.  Under  vacuum  conditions, 
the  cathode  was  shorted  to  the  anode  to  prevent  the  thruster  from 
eroding  while  it  was  fired  with  no  gas  pulse.  When  the  thruster 
was  fired,  magnetic  fields  similar  to  those  from  an  actual  fir¬ 
ing  were  created.  No  motion  was  indicated  on  the  dial  face, 
demonstrating  that  magnetic  perturbations  do  not  affect  impulse 
measurements  on  the  stand. 

A  similar  test  was  done  to  determine  air  drag  effects.  The 
calibration  stand  was  set  up  and  a  ball  was  dropped  several 
times,  as  is  normal  for  calibration.  The  chamber  was  then  pumped 
down  to  vacuum  and  the  ball  was  placed  on  the  holder  using  a 
remote  positioning  arm.  The  ball  was  dropped  the  same  number  of 
times  and  showed  the  same  deflection  as  that  for  open  air  test¬ 
ing.  This  was  as  expected,  given  the  small  wire  and  ball  sizes 
and  velocities  involved  (relative  to  system  mass) .  This  experi¬ 
mental  verification  confirmed  the  fact  that  open  air  calibration 
is  valid. 

Concerns  about  the  dial  indicator  interference  were  ad¬ 
dressed  by  setting  up  the  calibration  stand  and  a  measuring  rule, 
and  visually  measuring  the  thruster  displacement  with  and  without 
the  dial  indicator  attached.  The  thruster  displacements  were 
quite  small,  on  the  order  of  1/32  of  an  inch.  The  thruster 
motion  was  slightly  larger  with  the  indicator  removed,  but  how 
much  larger  could  not  be  determined.  It  was  not  more  than  3/64 
inch  total  displacement. 

Power  cable  interference  was  checked  (with  the  dial  indica¬ 
tor  attached)  in  a  similar  fashion  by  dropping  a  ball  with  and 
without  the  power  cable  attached.  The  difference  was  approxi¬ 
mately  one  division  on  the  dial  indicator  (approximately  l/lOOO 
inch) .  This  indicates  that  the  cable  provides  minimal  resistance 
to  thruster  motion.  With  the  small  velocities  involved,  any 
linear  damper  is  going  to  have  a  small  resistive  force.  Thus, 
the  cable  acts  as  a  very  light  damper,  with  almost  no  effect  on 
measured  thruster  displacement. 
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Table  4.  Measurement  Uncertainties 
(max  indicates  errors  may  be  smaiier, 
%  means  error  is  percentage  of  value) 


Item 

Units 

Range 

Uncertainty 

Ball  Horiz  Oist 

cm 

10-50 

+/-  0.5 

Dial  Displacement 

mils 

5-60 

+/-  1 .5  max 

Thruster  Mass 

kg 

58.127 

+5.  -1 

Plenum  Pressure 

psia 

20-30 

+/-  0.5 

Plenum  Temp 

deg  F 

80-100 

-►/-I 

Vacuum  Pressure 

Torr 

5E-S  -  5E-* 

+/-1% 

Thruster  Current 

kA 

0.8  -  30.0 

+/-  0.001 

Thruster  Voltage 

V 

50  -  750 

+/-1 

Gas  Pulse  Time 

msec 

106 

+/-  0.1 

Current  Pulse  Time 

msec 

1.5 -2.5 

+/-  0.01 

Local  Accel  Gravity 

m/sec 

9.79366824 

+/-1.8E-6 

Mass  Plow 

g/sec 

1.0 -2.0 

+/-1.2 

Mass  71ov  Errors 

Tho  mass  flow  is  a  crucial  da^um.  Mass  flow  mus't  be  known 
acctirately  and  precisely  to  get  realistic  performance  values.  The 
mass  flow  system  used  in  pulsed  MPD  testing  was  reasonably  pre¬ 
cise,  but  it  did  not  allow  an  accurate  mass  flow  measure.  The 
thermocouple  allows  only  an  approximate  temperature  measure, 
since  it  was  located  outside  of  the  chamber  instead  of  on  the 
plenum.  The  plenum-valve  system  made  it  very  difficult  to  attain 
the  same  pressure  for  a  test  set.  The  mass  flow  directly  depends 
on  the  plenum  pressure.  Difficulty  in  setting  the  exact  same 
pressure  was  not  an  error  in  itself,  but  to  speed  up  testing,  a 
variation  in  pressure  settings  was  accepted  by  the  test  team. 
This  led  to  various  mass  flows  being  taken  for  different  test 
cases.  The  "errors'*  arise  in  interpretation  of  the  data.  Which 
effscts  were  due  to  different  mass  flow  levels  and  which  were  due 
to  plasmadynamic  effects?  It  was  impossible  to  tell.  Finally,  the 
mass  flow  was  calculated  by  finding  the  pressure  change  before 
and  after  a  gas  pulse.  The  problem  was  that  the  pressure  change 
was  on  the  order  of  0.2  psi.  The  plenum  pressure  gauge  was  se¬ 
lected  for  expected  plenum  pressures  (order  30  psi) ,  so  that  the 
gauge  accuracy  was  on  the  order  of  0.5  psi.  Although  pressure 
changes  were  taken,  the  accuracy  of  those  measurements  is  uncer¬ 
tain,  as  is  the  absolute  accuracy  of  the  mass  flow. 

Well  after  the  test  matrix  was  conducted,  an  alternate  mass 
flow  calibration  was  devised  for  another  project  being  run  in  the 
pulsed  MPD  facility.  This  technique  consisted  of  flowing  propel¬ 
lant  into  another  plenum  in  place  of  the  thruster.  A  highly 
accurate  pressure  gauge  in  that  plenum  measured  the  change  in 
pressure  due  to  a  gas  pulse,  since  the  pressure  change  could  be 
measured  very  accurately,  using  the  ideal  gas  law  yielded  the 
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mass  flow  out  of  the  original  plenum.  Thus  a  correlation  between 
the  original  plenum  pressure  and  mass  flow  was  attained.  It 
turned  out  that  the  previous  mass  flow  calibration  curve  (ob¬ 
tained  using  the  old  technique)  was  correct.  Although  the  pres¬ 
sure  gauge  used  in  the  old  technique  was  used  below  its  stated 
accuracy,  the  detected  pressure  changes  were  good  enough  to  yield 
accurate  calibration  curves.  The  data  taken  with  the  old  tech¬ 
nique  is  probably  correct  in  mass  flow  magnitude,  but  since  the 
accurate  calibration  was  taken  long  after  the  data  was  taken,  the 
statements  made  in  the  previous  paragraph  must  hold. 

Voltagm  and  currant  Maasuremant  Errors 

Current  and  voltage  measurements  are  needed  for  power  effi¬ 
ciency  calculations  and  for  comparison  of  measured  thrust  to 
theoretical  thrust.  The  voltage  was  measured  by  an  attachment  to 
the  PFN  circuit.  Although  in  theory  the  location  should  give  the 
correct  thruster  voltage  (minus  a  minute  drop  through  the  power 
cables) ,  in  actuality  the  high  currents  involved  may  have  influ¬ 
enced  the  voltage  measurements.  All  the  testing  for  this  project 
was  done  with  voltage  probes  attached  to  the  PFN.  Recently,  the 
probes  were  reattached  at  the  power  feedthroughs  entering  the 
rear  of  the  chamber.  This  resolved  many  odd  current  and  voltage 
traces  previously  obtained  with  the  old  voltage  measurement. 
Unfortunately,  the  complete  test  matrix  had  already  been  taken 
with  the  bid  location,  resulting  in  erroneous  voltage  data  being 
recorded.  The  extent  of  the  uncertainty  is  not  clear,  but  it  may 
be  as  large  as  -f/-  15  per  cent. 

Another  concern  with  voltage  measurement  was  impedance 
matching  the  thruster  load  and  the  PFN.  This  concept  was  not 
considered  until  the  testing  was  over  with.  However,  impedance 
matching  is  an  important  step  in  the  MPD  testing.  If  the  thrust- 
er-PFN  system  is  not  properly  matched,  the  voltage  experiences  a 
ringing  or  underdamped  effect.  Although  this  effect  is  not  neces¬ 
sarily  deleterious  to  the  thruster,  it  does  invalidate  the  per¬ 
formance  results  for  comparison  with  other  researcher's  results. 
The  only  variable  changeable  in  the  MPD  facility  is  the  PFN 
inductance.  Although  it  is  not  very  convenient  to  change,  the  PFN 
inductance  can  be  made  to  match  almost  any  impedance.  The  imped¬ 
ance  is  a  function  of  PFN  current  and  capacitance,  both  of  which 
vary  with  the  level  of  PFN  voltage.  Voltage  variation  will  re¬ 
quire  matching  at  various  power  levels. 

The  current  was  measured  with  Pearson  coils.  These  coils  are 
actually  transformers  that  measure  electric  current  in  cables  run 
through  the  coils.  Although  the  coils  are  highly  accurate,  they 
do  have  a  saturation  limit.  That  theoretical  limit  is  at  a  PFN 
voltage  of  about  350  volts.  Since  testing  has  been  done  at  up  to 
750  volts,  it  is  clear  that  much  of  the  high  power  current  meas¬ 
urements  are  erroneous.  This  explains  the  odd  currents  and  volt¬ 
age-current  (V-I)  traces  at  the  high  power  levels.  Another  prob¬ 
lem  having  a  great  bearing  on  performance  measurement  is  the 
pulse  time  measurement.  This  time  pulse  is  used  to  calculate  the 
MPD  thrust.  Since  the  pulse  time  is  on  the  order  of  a  millisec- 
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ond,  slight  variations  in  the  measure  can  lead  to  large  varia¬ 
tions  in  calculated  thrust.  The  technique  used  to  estimate  the 
time  pulse  is  described  in  the  cold  gas  calibration  section.  With 
the  current  pulse  trace  tail  varying  in  position,  the  "equal 
areas  point"  has  varied  quite  a  bit.  The  variation  may  have  been 
due  to  Pearson  coil  saturation,  or  it  may  have  been  due  to  real 
effects.  The  end  result  has  been  comparatively  large  variations 
in  current  pul  a  time.  In  discussion  with  E.P.  researchers,  it 
was  agreed  that  the  best  estimate  of  the  pulse  width  is  at  the 
point  where  the  pulse  curve  has  dropped  to  95  per  cent  of  its 
maximum  height.  This  would  make  the  data  taken  at  PLEPL  more 
consistent. 


SUMMARY  OF  RESULTS 

Although  the  data  requires  more  verification,  it  does  pro¬ 
vide  trends  and  magnitudes  of  the  pulsed  MPD  thruster  perform¬ 
ance.  Unfortunately,  the  variations  are  difficult  to  attribute  to 
specific  factors.  As  previously  discussed  in  the  error  section, 
we  cannot  tell  whether  a  certain  curve  varies  because  of  thruster 
geometry,  differences  in  mass  flow,  improper  calculation  of  pulse 
time,  or  other  subtle  changes  in  test  configuration.  The  appendix 
contains  samples  of  voltage-current  traces  (the  most  accurate 
data  tad^en) ,  thrust  versus  power  (a  few  illustrative  examples  and 
a  summary  of  best  fit  curves) ,  specific  impulse  versus  power, 
efficiency  versus  power,  efficiency  versus  specific  impulse,  and 
thrust  versus  thruster  volume. 

Many  researchers  present  MPD  data  in  the  above  format,  with 
the  exception  that  some  researchers  compare  performance  versus 
current  rather  than  power.  The  PLEPL  results  compared  to  other 
researcher's  results  generally  indicate  a  worse  efficiency  for 
the  PLEPL  thruster.  Whether  this  is  due  to  low  accuracy  data,  due 
to  the  PLEPL  thruster  design,  or  due  to  the  test  configuration  is 
difficult  to  determine.  The  V-I  curves  did  compare  favorably  with 
those  reported  in  the  literature. 

The  final  item  in  the  list  of  graphs,  thrust  versus  thruster 
volume,  is  not  a  standard  method  of  presenting  data.  The  purpose 
of  this  graphical  representation  was  to  correlate  geometric 
magnitude  to  thruster  performance.  The  motivation  for  this  is 
that  MPD  thrust  is  produced  by  the  interaction  of  electric  and 
magnetic  fields.  Both  the  magnitude  and  orientation  of  these 
fields  determine  the  thrust  produced.  The  fields'  orientation  is 
directly  related  to  thruster  geometry.  To  a  first  order  approxi¬ 
mation,  the  thruster  volume  (open  space)  between  cathode  and 
anode  is  the  basic  geometric  variable  (assuming  that  the  anode 
and  cathode  shapes  do  not  vary,  but  the  electrode  lengths  vary) . 
The  three  types  of  data  points  (dots,  squares,  and  triangles)  are 
an  attempt  to  group  mass  flow  ranges  together. 

The  trends  indicate  that  thrust  is  inversely  proportional  to 
thruster  volume.  This  is  a  stronger  effect  as  power  levels  grow 
higher.  Assuming  the  data  is  accurate  enough  to  make  this  assess- 
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ment  (recall  the  discussion  in  the  error  section,  especially  on 
the  thrust  measurement  and  high  current  measurement) ,  we  can 
infer  that  for  smaller  volumes,  more  of  the  power  goes  into 
thrust  generation.  Perhaps  for  longer  electrode  lengths,  plasma- 
dynamic  instabilities  form  which  reduce  the  power  available  for 
propellant  acceleration.  This  result  needs  to  be  verified,  but 
the  implications  are  promising.  For  a  given  MPD  thruster  diame¬ 
ter,  the  smallest  length  (and  hence  weight)  may  be  the  best 
performing. 


LESSONS  LEARNED 

A  new  propellant  system  has  been  designed  for  the  MPD  facil¬ 
ity.  This  system  uses  a  fast  acting  piezo-electric  valve  and  a 
choked  orifice  to  provide  a  precise  and  accurate  flow  control 
system.  The  flow  through  a  choked  orifice  is  solely  a  function  of 
the  upstream  pressure.  Thus  instead  of  measuring  a  change  in 
pressure,  only  the  absolute  pressure  is  needed.  A  highly  accurate 
pressure  gauge  will  give  this  performance.  The  fast  acting  valve 
will  have  a  2  millisecond  open  and  close  time,  which  is  of  the 
pulse  time  order.  This  allows  a  very  rapid  settling  time  for  the 
mass  flow  and  requires  less  total  propellant  mass  ejected  during 
a  pulse.  It  also  allows  a  quicker  vacuum  chamber  pressure  recov¬ 
ery  time.  A  larger  plenum  with  an  attached  thermocouple  will 
allow  smaller  total  pressure  changes  and  a  direct  measure  of  the 
propellant  temperature. 

The  voltage  will  be  measured  solely  at  the  power  feed¬ 
throughs  leading  to  the  thruster.  Impedance  testing  will  be 
accomplished  to  determine  required  PFN  inductance,  and  future 
tests  will  involve  matched  impedance. 

The  data  will  only  be  taken  below  Pearson  coil  saturation 
limits.  Future  upgrades  will  include  purchase  of  Pearson  coils 
that  can  handle  currents  up  to  40  kAmp  with  little  saturation. 
Pulse  time  widths  will  be  taken  at  the  95  per  cent  drop  off 
point. 

Thruster  impulse  measurements  have  been  impeded  by  thruster 
mass.  The  thruster  was  designed  for  machining  ease  and  not  for 
optimal  experimental  configuration.  A  near  term  upgrade  will  be 
to  design  a  much  lighter  thruster  with  the  same  dimensions  of 
interest.  This  way,  the  thrust  pulses  will  cause  larger  displace¬ 
ments  that  can  be  more  easily  measured.  A  similar  pendulum  meas¬ 
urement  system  will  be  used,  but  a  redesign  will  eliminate  many 
of  the  old  system  uncertainties.  The  displacement  will  be  meas¬ 
ured  with  a  LVDT  or  the  acceleration  will  be  measured  with  an 
accelerometer,  thus  giving  a  direct  thrust  measurement.  The 
calibration  will  also  involve  a  similar  concept,  but  the  calibra¬ 
tion  frame  will  be  made  sturdier  and  an  improved  ball  release 
mechanism  will  be  used. 

Other  questions  remain  regarding  the  pulsed  MPD  facility.  A 
magnetic  field  characterization  would  be  helpful  in  understanding 
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possible  chamber  interactions  with  the  thruster.  Effects  of 
thruster  warm-up  on  thruster  performance  also  need  to  be  investi¬ 
gated.  And  finally,  some  thermocouple  measurements  on  heat  dis¬ 
tribution  over  the  thruster  might  also  be  helpful  in  overall 
understanding . 


FUTURE  EXPERIMENTS 

The  data  taken  thus  far  has  pointed  the  way  toward  an  im¬ 
proved  MPD  thruster  design.  However,  with  the  problems  encoun¬ 
tered  with  the  facility,  a  data  verification  must  be  run  with  the 
improved  facility.  The  verification  will  consist  of  a  detailed 
repetition  of  data  for  a  given  thruster  geometry  and  will  give 
a  measure  of  the  accuracy  of  the  previous  data.  In  addition,  an 
abbreviated  test  matrix  will  be  run,  repeating  several  of  the 
previous  electrode  combinations.  This  test  matrix  will  be  de¬ 
signed  using  the  Taguchi  method  to  get  maximum  performance  infor¬ 
mation  with  the  minimum  number  of  test  runs. 

Another  area  of  investigation  will  be  in  the  use  of  a  pre¬ 
ionization  plenum.  Part  of  the  motivation  for  this  investigation 
is  based  on  an  incident  which  occurred  during  the  initial  test 
matrix.  The  boron  nitride  backplate  came  loose  during  a  test 
firing.  The  thruster  motion  caused  the  plate  to  edge  forward. 
Since  the  plate  is  embedded  deep  in  the  thruster,  this  motion  was 
not  apparent  to  the  researchers.  But  the  sudden  jump  in  thruster 
performance  was  instantly  noticed  by  the  test  team.  What  had 
previously  been  ordinary  thrust  levels  and  current  pulse  shapes 
became  nearly  double  thrust  levels  and  pulse  curves  that  closely 
resembled  an  idealized  square  wave.  There  wasn't  time  to  further 
investigate  this,  so  the  PLEPL  will  investigate  a  pre-ionization 
plenum  at  a  later  date. 

A  possible  explanation  for  improved  performance  with  a  pre¬ 
ionization  plenum  is  that  the  PLEPL  MPD  uses  an  embedded  wire  to 
create  the  spark  that  breaks  down  the  propellant.  This  spark  is 
asymmetrically  distributed  through  the  MPD  cathode  channel,  as 
propellant  flows  and  expands  to  a  vacuum,  since  the  plasma  is 
sensitive  to  initial  conditions,  asymmetries  may  create  instabili¬ 
ties  that  reduce  MPD  performance.  The  pre-ionization  plenum  would 
provide  a  small  holding  area  for  better  propellant  ionization. 
This  ionized  propellant  would  uniformly  enter  the  main  MPD  chan¬ 
nel,  where  symmetric  fields  accelerate  the  propellant.  This 
conjecture  needs  to  be  verified  with  a  detailed  experiment. 


ERROR  ANALYSIS 

In  experiments  two  types  of  errors  propagate  to  cause  an 
uncertainty  in  the  final  results  (Ref.  12  and  13) :  systematic 
errors  and  random  errors.  Systematic  errors  are  those  that  occur 
due  to  faulty  calibration.  Random  errors  arise  from  the  inabili¬ 
ty  to  make  totally  accurate  measurements  from  one  observation  to 
another.  The  precision  of  an  experiment  depends  on  the  uncer- 
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tainty  due  to  the  random  errors.  This  section  discusses  the 
estimates  in  random  errors. 


Asstiming  a  Gaussian  distribution  in  the  experimental  values 
about  an  average  value,  the  uncertainty  in  determining  the  aver¬ 
age  value  is  proportional  to  the  standard  deviation  of  the  dis¬ 
tribution.  The  mean  value  of  the  distribution  and  the  "true” 
value  may  not  be  equal.  However,  as  can  be  shown  in  statistical 
theory,  there  is  a  95  percent  probability  that  the  difference 
between  the  "true”  value  and  the  mean  value  is  less  than  the 
uncertainty  given  by  the  standard  deviation.  This  defines  the 
uncertainty  in  the  results.  The  definition  does  not  give  a  value 
for  the  uncertainty,  but  it  does  allow  a  method  for  calculating 
the  uncertainty  of  the  results,  given  the  uncertainties  in  the 
parameters  that  are  used  in  the  mathematical  model. 

If  the  model  that  describes  the  system  has  at  least  two 
parameters , 

X  =  F(w,v, ...) 


and  the  most  probable  value  of  the  result,  X, 


is  given  by. 


X  =  F[u,v,,..) 


The  uncertainty  in  X  (ignoring  cross  terms)  ,  ,  is  given  by 

Equation  15 

or  for  Equation  9,  the  uncertainty  in  impulse  becomes 


Equation  16 


In  Equation  16,  the  O'"  's  are  the  values  listed  in  the  uncer¬ 
tainty  section  of  Tables  1  and  4.  The  values  of  the  partial 
derivatives  are  listed  in  Appendix  A.  with  the  parameter  uncer¬ 
tainties  determined,  a  computer  program  calculates  the  uncertain¬ 
ty  in  each  impulse  calculated,  given  the  measured  parameter 
values. 


CONCLUSION 

The  PLEPL  has  conducted  a  preliminary  experiment  in  the 
establishment  of  MPO  performance  for  different  electrode  length 
combinations.  Data  collected  has  given  us  an  understanding  of 
general  trends  in  thruster  performance,  but  limitations  in  exper¬ 
imental  capability  have  prevented  the  use  of  this  data  in  an 
absolute  sense.  The  most  valuable  outcome  of  this  project  has 
been  a  much  better  understanding  of  MPD  research  requirements  for 
the  PLEPL  staff.  The  pulsed  MPD  facility  can  now  be  upgraded  to 
gather  accurate  MPD  thruster  performance  data. 


24 


REFERENCES 


1.  Jahr,  R.G.,  Physics  of  Electric  Propulsion.  McGraw-Hill,  1968. 

2.  Kaplan,  David  I.;  Jahn,  R.G.;  Performance  Characteristics  of 
Geometr ica 1 Iv  Scaled  MPD  Thrusters .  MAE  1492,  Princeton  Universi¬ 
ty,  Feb  1982. 

3.  Mead,  F.B.;  Jahn,  R.G.;  Scaling  of  MPD  Thrusters .  Internation¬ 
al  Electric  Propulsion  Conference  Paper  79-2075,  Oct  30  -  Nov  1, 
1979,  Princeton,  New  Jersey. 

4.  Kaplan,  David  I.;  Jahn,  R.G.;  Performance  Characteristics  of 
Geometrically  Scaled  MPD  Thrusters .  MAE  1492,  Princeton  Universi¬ 
ty,  Feb  1982. 

5.  King,  D.Q.;  Jahn,  R.G. ;  Clark,  K.E.;  Maanetoplasmadvnamic 
Channel  Flow  for  Design  o£.  Coaxial  MPD  Thrusters .  MAE  1552, 
Princeton  University. 

6.  Phillips,  T.;  Bell,  K.;  100  Lbf  ACS  Engine  Evaluation  Program 
Facility  Checkout  Tests .  AFRPL-TR-79-01,  AFRPL,  Edwards  AFB,  Aug 
79. 

7.  Kaplan,  David  I.;  Jahn,  R.G.;  Performance  Characteristics  of 
Geometrically  Scaled  MPD  Thrusters .  MAE  1492,  Princeton  Universi¬ 
ty,  Feb  1982. 

8.  Mead,  F.B.;  Jahn,  R.G.;  Scaling  of  MPD  Thrusters .  Internation¬ 
al  Electric  Propulsion  Conference  Paper  79-2075,  Oct  30  -  Nov  1, 
1979,  Princeton,  New  Jersey. 

9.  Gilland,  J.H.;  Kelly,  A.J.;  Jahn,  R.G. ; MPD  Thruster  Scaling. 
International  Electric  Propulsion  conference  Paper  AIAA-87-0997 , 
May  11-13,  1987,  Colorado  Springs,  Colorado. 

10.  Sears,  F.W.;  Zemansky,  M.W.;  Young,  H.D.;  University  Physics. 
Addison-Wesley  Publishing  Co,  Reading,  Mass,  1983. 

11.  Kaplan,  David  I.;  Jahn,  R.G. ;  Performance  Characteristics  of 
Geometrically  Scaled  MPD  Thrusters.  MAE  1492,  Princeton  Universi¬ 
ty,  Feb  1982. 

12.  Bevington,  Philip  R. ,  Data  Reduction  and  Error  Analysis  for 
the  Physical  Sciences.  McGraw-Hill,  New  \ork,  1969. 

13.  Hogg,  R.V. ;  Tanis,  E.A.;  Probability  and  Statistical  Infer¬ 
ence.  MacMillen,  New  York,  1983. 


25 


APPENDIX  A 

Partial  Derivatives  for  Uncertainty  in  Mass  Flow 
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APPENDIX  B 


Partial  Derivatives  for  Uncertainty  in  Impulse 
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m  =  ball  mass 

M  =  thruster  mass 

L  =  string  length 

S  =  horizontal  distance 

k  =  coefficient  of  restitution 

Q  =  acceleration  of  gravity 
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APPENDIX  C 

Partial  Derivatives  for  Uncertainty  in  Thrust,  Specific  Impulse,  and  Efficiency 


Thrust,  T  1  =  Total  Impulse 
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APPENDIX  D 
Test  Matrix 


Note;  Lengths  are  in  centimeters.  Date  indicates  the  day  testing  began  on  an  anode-ca¬ 
thode  combination. 


Date 

Cathode 

Anode 

6  Sep  89 

3.2 

3.2 

8  Sep  89 

8.4 

8.4 

12  Sep  89 

3.4 

8.4 

15  Sep  89 

13.4 

13.4 

4  Oct  89 

8.4 

13.4 

11  Oct  89 

3.4 

13.4 

2  Nov  89 

18.4 

18.4 

9  Nov  89 

13.4 

18.4 

21  Nov  89 

8.4 

18.4 

27  Nov  89 

3.4 

18.4 

29  Nov  89 

3.4 

20.4 

1  Dec  89 

8.4 

20.4 

6  Dec  89 

13.4 

20.4 

8  Dec  89 

18.4 

20.4 

13  Dec  89 

20.4 

20.4 

1 1  Jan  90 

3.4 

3.4 

29 


DSA  601  DIGITIZING  SIGNAL  ANALYZER 
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Appendix  F 


Thrust  vs.  Power,  Specific  Impulse  vs.  Power,  Efficiency  vs. 
Power,  Specific  Impulse  vs.  Efficiency,  Thrust  vs  Thruster  Volume 
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Specific  Impulse  vs.  Power 
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Efficiency  vs.  Power 
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Efficiency  vs.  Power 
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Efficiency  vs.  Power 
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Efficiency  vs.  Power 
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Efficiency  vs.  Power 
11-21-89 


C:  8.4  A;  18.4 


99 


Efficiency  (%) 


Efficiency  vs.  Power 


11-27-89 


Power  (kWatts) 


C:  3.4  A:  18.4 


100 


0 


1 


1 

0  2000 


4-000  6000  8000  10000 

Power  (kWatts) 


12000 


Efficiency  vs.  Power 
12-8-89 


104 


Efficiency  (%) 


y 


Efficiency  vs.  Power 


0  2000  4000  6000  8000  10000  12000 

Power  (kWotts) 


105 


i 


Specific  Impulse  vs.  Efficiency 
9-6-89 


107 


Efficiency  (%) 


Specific  Impulse  vs.  Efficiency 


Specific  Impulse  (sec) 


108 


Lfficiency  {%) 


Specific  Impulse  vs.  Efficiency 
9-12-89 
30 

25 

20 

15 

10 

5 

0 

0  1000  2000  3000  ^000  5000  6000 


Specific  Impulse  (sec) 


109 


Efficiency  (%; 


Specific  Impulse  vs.  Efficiency 
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Specific  Impulse  vs.  Efficiency 
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Specific  Impulse  vs.  Efficiency 
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Specific  Impulse  vs.  Efficiency 
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Thruster  Internal  Volume  Geometry 


Note:  Shaded  region  denotes 
a  characteristic  volume,  not  an 
actual  physical  component.  The 
actual  thruster  is  a  copper  cylinder 
with  a  tungsten  rod. 
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Thrust  vs.  Volume 
For  2.0  Megawatts 


Volume  (cm''3) 

A  1.0-1.3g/s 
•  1.3 -1.7 
□  1.7 -2.0 
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